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Clump mass function at an early stage of molecular cloud 
evolution: II. Galactic cloud complexes 



o 

(N 

Oh 
< 

m 

(N 



< 

o 

6 



> 

O 
G\ 

in 

o 



* 



Todor V. Veltchev 1 ' 2 *, Sava Donkov 3 , and Ralf S. Klessen 2 

1 University of Sofia, Faculty of Physics, 5 James Bourchier Blvd., 1164 Sofia, Bulgaria 

2 Universitdt Heidelberg, Zentrum fur Astronomie, Institut fur Theoretische Astrophysik, Albert- liberie- Str. 2, 69120 Heidelberg, Germany 
3 Department of Applied Physics, Technical University, 8 Kliment Ohridski Blvd., 1000 Sofia, Bulgaria 



Accepted 2013, April 22 



ABSTRACT 

The statistical approach for derivation of the clump mass function (C1MF) developed 
by Donkov, Veltchcv & Klessen is put to observational test through comparison with 
mass distributions of clumps from molecular emission and dust continuum maps of 
Galactic cloud complexes, obtained by various authors. The results indicate gravita- 
tional boundedness of the dominant clump population, with or without taking into 
account the contribution of their thermal and magnetic energy. The C1MF can be pre- 
sented by combination of two power-law functions separated by a characteristic mass 
from about ten to hundreds solar masses. The slope of the intermediate-mass C1MF is 
shallow and nearly constant (—0.25 > Tim > —0.55) while the high-mass part is fitted 
by models that imply gravitationally unstable clumps and exhibit slopes in a broader 
range (—0.9 > Tim > —1-6), centered at the value of the stellar initial mass function 
(rHM < —1-3). 
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1 INTRODUCTION 

Dense clumps in molecular clouds (MCs) are typical sites of 
star formation as they are often associated with young stel- 
lar objects. Their origin can be sought in the early epoch 
of cloud evolution when supersonic turbulence creates sets 
of condensations in the cold, mainly molecular gas. Recent 
numerical simulations indicate their mean densities, sizes 
and masses vary in ranges 10 2 — 10 4 cm" 3 , 0.04 — 1 pc 
and 0.1 



Orion B JKramer et all Il998lh Taur us JOnishi et al.lll996t). 
Ophi uchus 1 Tachihara et alj 120001 ). Lupus ( Hara et al.l 
Il999l ) and many others. Some of these results on the 
physical parameters of clumps we r e inc luded in the 
statistical study of iTachihara et alj (J2002J) who consid- 
ered a sample of 9 Galactic star-forming regions. In 
the last decade, further intensive research of MCs by 
use of s ome high-den s ity t r acers like H 13 CO + (J = 
1 - 0) (lOnishi et all 120021: llkeda. Sunada fc Kitamural 



1LJMa Ti r S Ct i^ ^^XmX ' ^ B : llkeda. Kitamura fc Sunadal 1200911 as well o f dust 



120071 : iBaneriee etafl 120091 : IShettv et alj|2010l ). in consis 
tenc y with extensive observ ational data about MC clumps 
fe.g.. lBergin fc Tafallall2007l h Clump morphology is also di- 
verse: from filamentary to compact, quasi-spherical shapes 
(jHennebelle et alj I2008T ). Gravitational stability analysis 
shows that some clumps are subject to further contraction 
and collapse and eventually give birth to single stars or stel- 
lar clusters. || 

Numerous individual clumps were initially identified 
on maps of molecular line emissions which trace dif- 
ferent density regimes in MCs: 12 CO (n ~ 10 2 cm -3 ), 



Johnstone, Matthews fc Mitchell 



120051 [2006al: 

(e.g. lAlves. Lombardi 



continuum (IJohnstone et al. 2001 : Kerton et al. 2001 



_ 200J iReid fc Wilson] 
2010T ) and dust extinction 



Di Francesco et al.l. 

Ladal 120071 ) observations allowed 



CS (J = 1 - 0) and 
). Such surveys were 



13 CO and C ls O (n 

H 13 CO+ (J = 1 - 0) (n > 10 4 cm 

performed in nearby (at distanc es < 500 pc) Galacti c 

cloud complexes like Orion A (jTatematsu et alj Il993l ). 
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for more precise mapping of cloud structure. Some clumps 
originally found on emission maps were further decomposed 
and more compact (typical sizes < 0.2 pc), very dense 
(n > 10 5 cm~ 3 ) and probably collapsing clumps were 
delineated. Some authors call such objects 'dense cores'; 
hereafter, we label them simply cores. 

It is suggested t hat cores eventually form stars 
(|Bergin fc Tafallal 120071 ) and thus the study of their mass 
function (CMF) will enable a better understanding of the 
physical origin of the stellar initial mass function (IMF) and 
its possible variations. Indeed, numerous dust continuum 
and dust extinction observations demonstrate that the 
CMF resembles the IMF in its shape when fitted by a single 
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power- law JTesti fc Sargentl 1 19981: 1 Johnstone et alj|200ll). a 
combination of t wo pow er -law (IMotte. Andre fc Nerilll998l: 
Motte fc Andre! I200II; IJohnstone, Matthews fc Mitchelll 



20061 : 
tion ( 



Nutter fc W ard- Thompson 2007) or a l ognormal func- 
Stanke et al.l 120061 : lEnoch et all I2OO8I ; iKonvves et all 



ll998T ). Tau rus (lOnishi et al.l Il996f) . and of a sample of 
9 regions (jTachihara et al.l |2002i ); ii) from dust contin- 
uum surveys and CO map p ings of 2 kpc-away re gions: 
M17 jReid fc Wilson! l2006al: IStutzki fc Glisten! Il990l) and 



uOlOl l. It therefore has been proposed that the IMF is 
direct pro duct of the CMF and a unifo rm star formation 
efficiency (|Alves. Lombardi fc Ladall2007Q . 

Yet it is still unclear how the CMF originates from 
the mass distribution of the initially formed MC clumps. 
The clump mass function (C1MF), as derived from molec- 
ular line emission surveys of Galactic star-forming re- 
gions using low- and moderate-density tracers, varies in 
shape and slope and often differs substantially from the 
IMF. Most of the earlier studies of nearby complexes from 
CO mapping resulted in a single power-law C1MF of a 
rather shallow slope —0.3 > T > —0.8 without any char- 
acteristic mass Mrh JStutzki fc Giistenl Il990l: iBlit j Il993l : 



IWilliams. Blitz fc Stark|ll995l ; iHeithausen et al.lll998l) 
contrast to the Salpeter slope —1.3 ( Salpeterlll955T I of the 
high-mass IMF over M ch ~ O.5M . Moreover, the C1MF 
slope F was found to remain nearly constant for a large 
range of clump masses 1 0~ 3 ^ ui/Mq ^ 10 4 , in different 
MCs (|Kramer et alJll99SJ ). 

Later works allowed for a more detailed mapping of 
MCs. Emission from CO molecules was found to trace lower 
density cloud regions. Tracers like C ls O revealed structures 
with n ~ 10 4 cm -3 but essentially l arger ( I ~ 0.1 — 0.5 pc) 
than prest ellar cores (|Hara et al.l Il999l : iTachihara et al.l 
2000, 2002). Such compact clumps encompass about 10 C X of 
the cloud mass. It was demonstrated that the derived CIMFs 
could be represented by a combination of 2 power-law func- 
tions as M c h ~ 15 — 20 Mq separates intermediate-mass 
from high-mass C1MF, with slopes —0.3 > Tim > —0.7 and 
Fhm < —1.5, respectively. This general picture was con- 
firmed from submillimeter continuum studies of MC com- 
plexes at kpc-scale distances. Using 450 fim and 850 /im 
SCUBA maps, Reid & Wilson derived 2 power-law CIMFs 
for two massive star-forming regions. For NGC 7538, they 
obtained Af c h ~ IOOMq and Fim ~ while the vari- 
ations of the high-mass slope seem to be more sensitive 
to the wavele ngth and close to th e Salpeter value: — 1 > 
Fhm > -1.6 (JReid fc Wilson! 120051 ). For M 17, the charac- 
teristic mass turned out to be an order of magnitude less 
(M c h ~ 8 — 20 Mq), the intermediate mass slopes are posi- 
tive (Fim ~ 0.5 — 1) and the high-mass one s are about the 
value for CO clumps -0.5 > r HM > -0.9 <|Reid fc Wilson! 
l2006al ). Extending further their study to 11 low- and high- 
mass star-forming regions, those authors conclude that a 
single power-law C1MF is clearly ruled out and argue for a 
two power-law function with —0.2 > Tim > —0.7 and Thm 
about the Salpete r slope, regardless of th e diversity of char- 
acteristic masses (jReid fc Wilsonll2006bl ). 

The predictions of our statistical model 
IjDonkov. Veltchev fc Klessen! 120121 . hereafter, Paper I) 
are generally consistent with a two power-law C1MF 
with shallow intermediate-mass slope Fim and steeper 
high-mass slope Thm, allowing for a variety of Af c h within 
two orders of magnitude. In this work they are com- 
pared with results on the observational C1MF: i) from 
molecular -line studies of three ne arby star-fo rming regions: 
Orion A (jTatematsu et al.lll993T ). Orion B (jKramer et al.l 



Rose tte ( Di Francesco et all bold : IWilliams, Blitz fc Starkl 
Il995f l. We define clumps as condensations formed through 
a turbulent cascade during the early MC evolution. The 
ensemble of clumps generated at given spatial scale L 
obeys a power-law relationship between clump masses and 
densities n oc m x where the exponent x = x(L) is calculated 
considering equipartition relations between various forms of 
energy: gravit ational, turbulent (kinetic), therma l (internal) 
and magnetic IjDonkov. Veltchev fc Klessenll201ll . hereafter, 
DVK11). The C1MF is derived as a superposition of clump 
mass distributions over a range of scales, taking into 
account the fractal structure of the cloud (Paper I). 

The physical basis of the model, its free parameters and 
the construction of the C1MF are recalled in Section [2] The 
clump mass distribution derived from modeled structure of 
each considered individual cloud as well as different distri- 
butions that fit the composite observational C1MF are pre- 
sented and commented in Section [3] Section [4] contains a 
discussion on the applicability and the restrictions of our 
approach to predict the observational C1MF and envisions 
its possible extensions. Our conclusions are summarized in 
Section [5] 



2 STRUCTURE AND CLUMP MASS 
FUNCTION OF INDIVIDUAL CLOUD 

2.1 Physical framework of the model 

Our statistical model for clump description is presented in 
details in DVK11 (Sect. 2 and 3) and Paper I (Sect. 2). Its 
main assumptions can be summarized as follows: 

(i) Scaling laws within the turbulent cloud: We consider 
fully developed supersonic turbulence that implies homo- 
geneous and isotropic stochastic medium with a fractal 
structure and well-defined scaling laws of turbulent veloc- 
ity, mean density and mean magnetic field. Turbulent flows 
create density structures at any scale in the inertial range 
L max > L > I/min through a casca de possibly driven by 
the v ery process of cloud formation IjKlessen fc Hennebelld 
I2OIO1 ). We estimate the upper limit I/ max = 20 pc adopting 
a typical size of giant MC ~ 50 pc as an injection scale and 
taking in view that the largest scale o f the turbulence in- 
ertial range is abou t a factor of 3 less (jKritsuk et al.l 120071 : 
IPadoan et al.ll2006l ). The lower limit L m i n is imposed above 
the actual end of the inertial range from the construction of 
our model: the sizes I of all clumps generated at given scale L 
must be within the inertial range and are typically an order 
of magnitude less than L. Various observations show that 
the inertial range spans at least 3 orders of magnitude, i.e. 
the size of the smallest generated object l m i n must be about 
0.02 pc. We consider mainly molecular and isothermal gas 
with temperatures T = 10 — 20 K. Requiring supersonic 
medium at all fractal scales and by use of a typical velocity 
scaling (see equation [T] below), one obtains l m i n > 0.05 pc 
which yields a scale of its generation L m i n ~ 0.5 pc. 

Turbulent velocity dispersion u and mean mass density (p) 
are assumed to scale according to "Larson's first and second 
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laws" (|Larsonlll98ll ) whereas the scaling relation of the mean 
magnetic field B is obtained from its observational ly verified 
relati on to the mean mass density (B oc (p) 1 ' 2 : ICrutcherl 
1 19991 ): 



U=1.1L P [kms -1 ] 
(p) = 13.6 x 1(T 21 L° 



[gem 



B = 50 L" 



[pG] 



(1) 
(2) 
(3) 



A fixed velocity scaling index /3 is chosen. While observa- 
tional and numerical studies indicate a broad range 0.2 < 
P < 0.65, the results in this Paper narrow it to 0.2 < 
/3 ^ 0.47 ('soft' velocity scaling). The density scaling in- 
dex a = a(L) is derived self-similarly from the assumption 
of mass-density relationship for clumps generated at a given 
scale L (see below). 

(ii) Lognormal clump density distribution: Such volumet- 
ric distribution of mass density p is testified from numerous 
numerical simulations of supersonic turbulence and is de- 
scribed through a standard lognormal probability density 
function (pdf): 



p(s) ds 



1 



V2~7 



exp 



ds 



(4) 



where s is the log density and (s max , a) are the distribution 
peak and the standard deviation, respectively. In our model, 
this density statistics is used as the basis of clump statistics 
as follows. Let N^ ot is the total number of pixels in a cloud 
map and Np is the number of those with density p. Then 






1 



V27] 



exp 



Since turbulence is assumed to be homogeneous and 
isotropic and gravity is a central force, one can consider 
the generated clumps as homogeneous spheres, character- 
ized solely by their size (diameter) I. Of course, the real 
fragments in a turbulent MC with high Mach numbers (^ 3) 
delineated by isodensity contours [p — Ap, p + Ap] can be 
with complex shapes and even not necessarily connected re- 
gions. The essence of our statistical approach is to describe 
them by an ensemble of Np oc (Np /N£ ot ) spherical homo- 
geneous clumps with density p. Thus the 'average clump 
ensemble' generated at given spatial scale L follows a log- 
normal density distribution like the pdf at that scale. 

The parameters of the density distribution: 



a 2 = \n(l + b 2 M 2 ) , s m ax = -y 



(5) 



depend on the spatial scale through the sonic Mach number 
.M = u(L)/cs where c s is the sound speed. The turbulence 
forcing parameter b spans values between 0.33, for purely 
solcnoidal forci ng, and 1.0 0, for purely compressive forcing 
(see lFederrath et all (|2010l )). The modeled CIMFs presented 
in this Paper favor mainly solenoidal forcing, mixed in a 
natural way with compressive modes: 6 ^ 0.46. A comment 
on that is included in the Discussion. 

(iii) Clump mass-density-size relationship: Masses, densi- 
ties and sizes of clumps in an 'average ensemble' are assumed 
to obey the statistical relationships: 



In 



\m J 



h 'sJ-T 



3.C 



Wo 



with a choice of normalization units: 



po 

h 

Poll 
m 



(p) 
kL 



2 1 

exp | a x 



(6) 
(7) 



(8) 
(9) 

(10) 



The dimensionless parameter k accounts for the precision 
of the clump size 'measurement' and can be interpreted as 
mapping resolution of the scale of clump generation. It is 
appropriate to set it as a small constant of order of several 
percent since I is typically an order of magnitude less than 
L. 

According to our turbulent scenario of clump formation, 
the clump density-size relationship (equation]?]) is taken as a 
self-similar extension of the density scaling law (equation [5| 
and hence: 

a = JE_. (ii) 

1-x v ' 

(iv) Equipartitions between clump energies: A general 
type of equipartition relation between gravitational W and 
kinetic £kin energy per unit volume with some coefficient 
/ g k of proportionality 

|W|~/gk-EW, 

is expected to hold for structures shaped by turbulence in 
which gravity gradually takes over; e.g. in regions where 
turbulence decays locally or where it accumulates material 
reaching a state of local gravitational instability. A fidu- 
cial range 1 ^ / g k ^ 4 is testified from simulations of 
the early stag e of the clump evolution, befo re stars have 
been formed (|Vazauez-Semadeni et al.l 120071 ). The analy- 
sis of MHD simulations of cloud formation at Galactic 
scale (length unit ~1 kpc: iPassot et al.lll995l ) showed that 
the considered equipartition can include contributions of 
thermal (internal) _E t h and magnetic J? maE ene rgy as well 
(jBallesteros-Paredes fc Vazquez-Semadenilll995l ). 

In our model, we assume such equipartition relations to 
hold for the 'average clump ensemble' and use them to derive 
the mass-density scaling index x(L) at each scale: 

• wkin2 or wkin4, equipartition of the gravitational vs. 
kinetic energy: 



\W\ ~ 2E kin 
\W\ ~ 4S kin 



(12) 
(13) 



• wkin2mag, equipartition of the gravitational vs. ki- 
netic and magnetic energy: 



\W\~2Eun + E„ 



(14) 



• wkin2th2, equipartition of the gravitational vs. ki- 
netic and thermal energy: 



\W\ ~ 2E Mn + 2E t: 



(15) 



The corresponding equations for the typical member of 
the 'average clump ensemble' which are used to derive 
x(L) are listed in Paper I (Appendix A and B). 
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2.2 Free parameters of the model and cloud 
structure from extinction maps 

In our model, a MC is considered as a hierarchical set of 
spatial scales and its structure is described through the so- 
lutions x(L) obtained for a chosen equipartition relation. 
They depend on four free parameters: velocity-scaling in- 
dex p, turbulent forcing parameter b, mapping resolution 
k (equation [9]) and temperature T. The parameters' ranges 
of variation that yield plausible solutions can be restricted 
b y use of observational studies of MC structure. The work 
of lLombardi. Alves fc Ladal (|2010l . hereafter, LAL10), based 
on dust extinction maps of Galactic cloud complexes, corre- 
sponds best to our approach. These authors define an 'effec- 
tive radius R B = \/ S/-k of a subregion or a set of subregions 
with total area S which we interpret as an observational 
counterpart of our notion of 'spatial scale', i.e. L = 2R B . 
Then the MC structure is described by the relationship be- 
tween the effective radius and the mass within the total area 
S: M B oc R] (see Fig. 2 in LAL10). On the other hand, the 
total mass contained within a scale L depends, in our ap- 
proach, on the index x(L) due to the self-similarity assump- 
tion (cf. equations [7] and EJ: M(L) = (p)L 3 oc L 3 ^ 1 -^. 
Hence, the mass-scale power-law index 7 = 3/(1 — x) or: 



(16) 



Ml 



1 1- 



Comparison between x(2R s ) obtained in that way from 
the work of LAL10 and our model predictions x(L) for sets 
(j3, b, k, T) is described in Paper I (see Fig. 1 there). The 
best fits x(L) for the cloud complexes Taurus, Orion A and 
Orion B studied in this Paper are plotted in the left panels 
of Fig. [T] [2] and [3 respectively. We comment them in Section 
® 



2.3 Derivation of the C1MF 

A parameter set (/3, b, n, T) corresponding to the obtained 
best fit x(L) of the LAL10 data for a chosen MC yields: i) a 
lognormal clump mass distribution pl (m) ; and ii) a measure 
of the total number of clumps N to t(L) at each scale (DVK11, 
Eq. 17): 



, r , r . 1 / 2 (z-l)(l-2a;) 
N tot (L) = -, exp I a 2 x i A_ >- 



(17) 



The C1MF of an individual cloud is derived as a super- 
position of the lognormal clump mass distributions. Each 
of them is discretized through sets of weights {NL(mj) — 
PL(mj)N t ot(L), j = l,...,n} where the mass range is cen- 
tred at the peak of the distribution and the mass limits are 
determined from the requirement: 



N t , 



t (L) = X> L (m,) 



(18) 



Also, to construct the C1MF correctly, one must re- 
quire inas^_conservation throughout the MC fractal struc- 
ture (|Elmegreenl 1 1997T ) . The whole cloud with size equal to 
the upper limit of the inertial range i m ax = 20 pc is to con- 
tain Ml substructures of mass M(L) at given scale L. As 
shown in Paper I (Sect. 3.2): 



£!-*(£) 



(19) 



Eventually, one obtains for the C1MF value in a selected 
mass bin m! — Am' ^ m ^ m' + Am': 



FciMp(m') = ]T Ml ]T N L (m) 



(20) 



If only gravitationally unstable clumps are considered, 
one should take also into account their contraction in 
timescales given by the free-fall time rg oc p -1 . Their 
time-weighted mass function can be derived by permitting 
constant replenishment of the clump population and intro- 
ducing a weighting factor of each scale of clump generation: 
w(L,x(L)) oc TgT 1 oc (p) 1/2 oc £3^/2(1-*) ^ equations [2] 
and llip . Then the time- weighted C1MF is obtained as the 
number of substructures Ml in equation [20] is modified by 
factor w(L, x(L))/w(L max , x(L max )). 



3 RESULTS 

We put to test the C1MF, derived in Paper I, using molecu- 
lar line observations of nearby Galactic MC complexes and 
dust continuum studies of further 3 regions at kpc-scale dis- 
tances, mapped also in CO lines (Table [l]). Subsamples of 
the original clump data are considered taking into account 
the lower mass limit of confidence in our models and/or the 
observational completeness limit; the corresponding size and 
mass ranges are specified in columns 3 and 4. 



3.1 Molecular line studies 

Three individual nearby star-forming regions of comparable 
size are selected: i) with clumps deli neated through sim - 
ilar density tr acers: 13 CO (Orio n B; iKramer et al.lll998T ). 
C 1 8 Q (Taurus;IOnishi et allll996h and CS (J = 1-0) (Orion 
A: lTatematsu et al.lll993T ). consistent with the typical clump 
densities in our approach n ~ 10 3 — 10 4 cm -3 ; ii) with known 
general structure from dust extinction maps of LAL10. The 
spatial extent of the studied regions in the chosen com- 
plexes falls within the adopted inertial range of turbulence 
(Sect. [2~T1 (i)). Dust extinction data reveal some diversity 
of structure in terms of mass-size relation of regions de- 
lineated by isodensity contours: Taurus exhibits monoton- 
ically 'steep' structure in terms of x(L) while Orion A and 
Orion B have 'shallow' internal regions and 'steep' exter- 
nal, less dense regions (Fig. IU3I left panels; cf. Fig. 2 in 
LAL10). Taurus is a typical low-m ass star-forming region 
(jKenvon. Gomez fc Whitney! 1 20081 ) . while numerous high- 
mass stars have been formed in Orion A and B in the recent 
past. As an additional test of the predictive power of the 
models for a variety of environments, we perform compar- 
isons with observational C1M F for a sample of 9 MC com- 
plexes (jTachihara et al.ll2002f ) wherein star-forming, cluster- 
forming and starless clumps have been detected. 

An equipartition relation and a parameter set were 
sought that yield the best fit x(L) of the cloud structure 
as traced by the LAL10 data (Fig. [QO left panels). The 
modeled C1MF derived from them is compared with the ob- 
servational one from molecular line maps (Fig. IU3I right 
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wkirvt, p=0.33, b=0.40, T=10 K 
k=0.075 - 
k=0.085 




100 



t^=H4 



-0.27 + 0.17 
-0.46 ± 0.09 



100 



L [pc] 



m I M Syn I 



Figure 1. Structure and C1MF of the Taurus complex. Left: Best fits x(L) of the data from LAL10 (dots with error bars) are plotted 
with lines and their characteristics (equipartition relation, model parameter values) arc specified. Right: Clump mass functions, derived 
from the same models (lines), are juxtaposed with the observational C1MF from molecular line maps (open symbols with error bars). 
The range of scales of clump generation (left) corresponding to the observational clump mass range, restricted by the lower limit of 
confidence of the model (right) arc shown with shaded areas. A weighted power-law fit of the observational C1MF (blue line) is drawn. 
Its slope (obs) and the one of the modeled intermediate-mass C1MF (mod) are specified. 



Table 1. Observational clump samples used for comparison with the predictions of our model. Regions that are studied both by use 
of molecular line emissions and dust continuum data are put in bold. The slope estimates are re-derived by us, adopting mass limits 
corresponding to the applicability of our model and to the data completeness range. Notation: D = approximate distance to the region, 
IM = Intermediate-mass, HM = High-mass, -y obs = Slope [d log m/d log r] on the clump mass-size diagram, M c ^ = characteristic mass, 
CF = Clump-finding (method), pHM = contour at Half-Maximum around intensity peaks, pPV = Position- Velocity diagrams passing 
through intensity peaks; GB = gravitationally bound, UB = unbound 



SF region 



D 

kpc ; 



Ref. 



Sizes 

[PC] 



Masses 
[Mo] 



"Yobs 



[M Q ] 



C1MF slope 
IM HM 



CF method 



Note 



Molecular line studies 



Taurus 
Orion A 
Orion B 

MC sample 

M17 
Rosette 



0.14 
0.45 
0.45 
^ 0.18 
^ 0.18 
1.6 
1.6 



Dust continuum studies 
M17 1.6 7 

M 17 + NGC 7538 7,8 

Rosette 1.6 9 



0.07 - 0.45 
0.06-0.35 
0.06-0.70 
0.08-0.45 
0.08-0.45 
0.03 - 0.60 
0.45-4.5 



0.02 - 0.70 
0.06 - 0.60 
0.20 - 0.70 



3-80 
30 - 1000 

3-200 
3-60 

3- 100 
6 - 2300 
30 - 2000 



1-600 
3 - 1000 
1.7-200 



1.95 
1.14 
1.95 
2.09 
2.09 
1.77 
2.40 



2.13 
2.05 
2.22 



^ 100 
-300 
~ 13 
~ 10 
~ 13 
> 200 
>800 



- 100? 

- 100 
>9 



-0.27 
-0.14 
-0.17 
-0.25 
-0.44 
-0.34 
-0.28 



-0.35 
-0.25 
-0.57 



-1.82 
-1.28 
-1.45 
-1.77 
-1.00 
-1.13 



-1.82 
-1.04 



pHM 

pPV 

Gaussclumps 

various 

various 

Gaussclumps 

Clumpfind 



Clumpfind 
Clumpfind 
Getsources 



mostly GB 
~90% starless, GB? 

UB + GB 
mainly starless, GB? 
mainly starless, GB? 

UB + GB? 
starless, mostly UB? 



mainly starless 
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panels). The range of spatial scales considered in the fit- 
ting procedure (shaded areas, left panels) corresponds to the 
mass range of the referred observational C1MF (shaded ar- 
eas, right panels), additionally restricted by the lower clump 
mass limit of confidence in the model. The scales are treated 
as those where the observed clumps were generated and 
hence their range is obtained through the clump mass-scale 
diagram of the tested model (Paper I, Sect. 3.1). 



3.1.1 Taurus 

Fig. [T] (left) demonstrates that the equipartition choice 
wkin4 yields a very good fit of the inner parts of the 
complex. That seems physically consistent to us since this 
case describes structures which have been shaped under 
significantly influence of self-gravity. The typical dense 
clumps generated in them have mass distributions in a rela- 
tively good agreem ent with the o bserv ational C1MF (Fig. [l] 
right), obtained bv lOnishi et al.l (J1996T ) ■ The best-fit param- 
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Figure 2. Structure and C1MF of the Orion A complex. The designations are the same like in Fig. [T] slopes of both intermediate-mass 
and high-mass C1MF are specified. The modeled non-time-weighted (lines) and time-weighted CIMFs (hatched area) are shown. The 
completeness limit of the observational data is adopted as a lower mass limit of confidence. 



eters point to a mainly solenoidal turbulent forcing (see 
iFederrath. Klessen fc Schmidtl 120081 ) and a 'soft' velocity 
scaling, typical for incompressible turbulence. On the other 
hand, the equipartition relation wkin4 is obviously not ap- 
plicable to the less dense extensive parts of the cloud. We 
stress, however, that in our model scales L > 3 pc do not 
produce clumps with masses within the considered observa- 
tional range. 

It should be p ointed out that the C1MF obtained by 
lOnishi et al.l |l99y) differs significantly from CIMFs in other 
clouds as derived from molecular line emission data. Using 
essentially the lowest mass bin (~ 3 M@), these authors 
argue for a flat mass function^] in Taurus, although they 
remain open for an alternative interpretation based on an 
"unbiased survey... in a wide mass range". We believe that 
indeed their data cover the mass range around the turn- 
over of the C1MF while its high-mass part (if existing) falls 
beyond their scope - note that characteristic masses M c h of 
hundreds Mq are derived in various MC complexes (Table 1, 
column 7). The power-law slope of the modeled C1MF over 
the lower mass limit of confidence in our model is slightly 
shallowe r than the typical one of the mass distribution of CO 
clumps (|Blitdll993l : iKramer et "all 1 19981 ). It would steepen 
(r ~ —0.7) if the observational data contained clumps with 



A shallow mass spectrum, in their terms. 



masses up to the expected Mch ~ 200 Mg 
CIMFs from wkin4 (Fig. 3 in Paper I). 



3.1.2 Orion A 



cf. the modeled 



We obtained good fits of the Orion A structure in the case 
wkin2mag at intermediate to large spatial scales (Fig. [21 left). 
The corresponding CIM Fs generally agree with the result of 
iTatematsu et al.l (|1993f >. Enhancement of the model tem- 
perature, conrp_atlble_wMi_tlie_de£^ed_T ~ 20 K in the 
complex (jlkeda. Sunada fc Kitamurall2007l ). leads to larger, 
more plausible values of b and j3 and yields a wider range of 
scales of clump generation as the cloud structure is better 
fitted through the curve x(L) (Fig. [2] bottom). 

Depending on the model temperature, the case 
wkin2mag in our approach generates mostly/only gravita- 
tionall y unstable clump s (Pap er I, Fig. 3). On the other 
hand, ITatematsu et al.l (|1993T ) claim that all clumps in 
their sample are close to virial equilibrium. Therefore a 
time- weighted C1MF (cf. Sect. [273]l is more appropriate for 
comparison with the referred observational study. Time- 
weighting steepens the high-mass slope of the model from 
Thm 1 ~ — 1> typical for fractal clouds, to ~ —1.5, i.e. steeper 
than that of the stellar IMF (F ~ -1.3; ISalpeterl H955T ) 
(Fig. [2] bottom right). This leads to agreement (within the 
data uncertainties) with the observational C1MF for clump 
masses > M c h ~ 300 Mq while for lower masses the devia- 
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Figure 3. Structure and C1MF of the Orion B complex. The designations are the same like in Fig. [2] 



tion is drastic. In contrast, the non-weighted modeled C1MF 
is generally consistent with the observational C1MF below 
M c h and for T = 20 K. We note the high completeness limit 
(~ 30 M©) of the data which might artificially constrain 
the intermediate-mass range. 



3.1.3 Orion B 

The ref erred observational C 1MF in this region is part of the 
work of lKramer et al.l (|l998T ) who analyzed CO data sets for 
8 Galactic MCs. They argued for a single power-law CIMFs 
in all cases, with a universal mean slope F ~ —0.7, over a 
wide range of clump masses. However, a closer look at their 
result for Orion B reveals that this MC may be an exception. 
The conclusion of lKramer et al.l (|1998l ) about a single power- 
law with shallow slope is based mostly on the single mass bin 
> 200 Mo (Fig-El r ight). Inspection of the statistically rich 
data in the range 2 — 80 Mq indicate rather a combination 
of two power laws: a shallow intermediate-mass slope below 
15 Mq and a significantly steeper high-mass slope. Adopting 
the latter value as M c h and the completeness data limit of 
the authors M ~ 1 Mq, one obtains Fj^f similar to that in 
Orion A and a Salpeter-like Fh'm- 

Modeling of the Orion B structure as traced by the 
LAL10 data yields good fits in the cases wkin2th2 and 
wkin2 (Fig. [3J left). For such choices of equipartition and 
temperature T — 10 K, the critical mass over which most 
of the generated clumps are gravitational unstable is about 
the adopted characteristic mass from the used observational 
data. Therefore we apply time-weighting for the high-mass 
regime and obtain T^m 1 close to the Salpeter value. The 



modeled CIMFs fit remarkably well the high-mass observa- 
tional ones, excluding the sole bin at the upper limit of the 
distribution (Fig. [3J right). Non- weighted CIMFs obviously 
fail to reproduce the observational one altho ugh the high- 
mass slope ~ —0.7 is close to the estimate of lKramer et al.l 
l|1998r ). The model is not applicable to fit the intermediate- 
mass C1MF because of the high lower mass limit of confi- 
dence. 



3.1.4 Sample of MC complexes 

The presented best-fit descriptions x(L) of MC structure 
in Taurus, Orion A and Orion B and their corresponding 
CIMFs with different characteristic masses M c h and HM 
slopes (cf. Table [1} are derived by use of different equipar- 
tition relations that may reflect a variety of physical con- 
ditions. What model case would describe appropriately a 
statistical C1MF obtained from a large set of clump data in 
differe nt MCs? The sample in the study of iTachihara et al.l 
(j 20021 ) could give a clue to the answer - it encompasses 
nearby cloud complexes of diverse star-forming type as the 
vast majority of clumps are associated with no or a small 
number young stellar objects and have masses below sev- 
eral dozens Mq. Clumps associated with young stellar clus- 
ters were found in two out of nine sampled clouds and their 
fraction is only 5 % of the total number of clumps. All ob- 
jects with masses ^ 115 Mq are from this group. One may 
take the latter value as an upper mass limit when fitting the 
C1MF, in order to avoid extreme environments and to secure 
sufficient statistics. The two power-law shape of the C1MF 
is evident whereas the characteristic mass is ambiguous. We 
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Figure 4. Modeling of the observational C1MF from 9 Galac- 
tic star-forming regions by use of the equipartition relations 
wkin2 and wkin2th2 for k = 0.065 and adopting characteristic 
mass ~ 10 Mq (top) and ~ 13 Mq (bottom). See text for the 
considered mass ranges. Other designations are the same like in 
Figs. HE] 



consider two choices which both allow fo r very good fits: 
i) M c h ~ 10 Mq, adopting the estimate of iTachihara et al.l 
(J2002f ) and decreasing the upper mass limit further by one 
bin; and ii) M ch ~ 13 M (Fig. g]). The lower mass limit 
in both cases is the lower limit of confidence of the tested 
models. 

The choice i) yields excellent agreement of the model 
case wkin2th2 with the observational C1MF: time-weighted 
model for the HM part and non-time-weighted for the IM 
part (Fig. 2| top). The IM slope is very similar to the ones 
found in the considered individual MC complexes while 
the HM slope is about the Salpeter value, like in Orion 
B (Sect. I3.1.3[) . The discrepancy with the data for m > 
60 Mq could be attributed to incompleteness of the sam- 
ple in this mass range. A Salpeter-like slope is to be ex- 
pected for mass distributions of (nearly) virialized clumps 
assuming their formation at a constant rate. The effect will 
be similar if the C1MF is derived from a statistically sig- 
nificant sample containing clumps in complexes at different 
evolutionary stages of clump formation. Note also that the 
variation of the velocity scaling index /3 in the fitting models 
is co nsistent with most observ ational and numerical works 
(e.g. IPadoan et alll2006l , |200Qt ). The values of the b param- 



eter indicate preliminary solenoidal turbulent forcing with 
small contributions of the compressive mode. 

Considering the choice ii), an upper mass limit ~ 
115 Mq is adopted as mentioned above. The model case 
wkin2 provides an excellent fit of the intermediate-mass 
C1MF and a problematic one - for the high-mass C1MF. 
Again, the completeness of the clump mass bin 60 ^ m ^ 
100 Mq is an open issue which is crucial for the correct 
slope estimation. On the other hand, the obtained low ve- 
locity scaling index f3 is consistent with results for Taurus, 
Orion A and Orion B. 

Modeling of the observational C1MF from the equipar- 
tition cases wkin2 and wkin2th2 suggests that most of the 
clumps generated in MC complexes of diverse star-forming 
activity and with mean density about the typical local den- 
sity (cf. Eq. 7 in Paper I) are in a state close to virial equi- 
librium. 



3.1.5 Clump mass-size diagrams 

An additional test for adequacy of the proposed statisti- 
cal modeling of the C1MF is to construct clump mass-size 
diagrams. The comparison between the used observational 
data and the best-fit models of the observational CIMFs 
is displayed in Fig. [5] There are two difficulties with such 
analysis. First, the 'average clump ensembles' in our model 
occupy narrow strips on m— r diagram^ whereas the disper- 
sion of the observational samples is huge. Second, mass and 
(especially) size determinations from molecular-line maps 
depend essentially on the clump-finding technique which is 
different in each considered reference work (Table [l] col- 
umn 10). A meaningful criterion is to compare the slope 
7° bs = (dlogm/dlogr) on the observational m — r diagram 
with the predicted slopes 7 mod as the latter vary with the 
chosen equipartition case and the velocity scaling index p 
(cf. Paper I, Fig. CI). 

The predicted slopes from the best-fit C1MF model 
are consistent with the ones from the referred molecular- 
line studies within la limit for Taurus and Orion A and 
within about 2<t limit for Orion B and for the sample of 
ITachihara et al.l (|2002l ). In two cases, the 'average clump 
ensemble' strips cross the region of the observational data 
(Fig. [5] bottom). Taking into account the good agreement 
between observational and modeled CIMFs (Figs. [TM)| . the 
shift between the two sets on the m — r diagrams - with a 
factor of 2 toward larger (Orion A) or smaller sizes, - is to be 
interpreted with the variety of the clump size definition. The 
issue will be commented further in the next Section wherein 
we compare our model with CIMFs from dust continuum 
studies and in the Discussion. 



3.2 Dust continuum studies 

Sensitive submillimeter continuum maps are appropriate for 
probing the MC structure on small scales due to the op- 
tically thin dust emission. In the last decade, the achieved 
angular resolutions of < 20" from ground-based observations 
(e.g. SCUBA, Bolocam) and space missions like Herschel 



2 Clump sizes from our model are converted to radii r = 1/2 to 
be compared with the radii of observed clumps. 
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Figure 5. Clump mass-size diagrams from the referred molecular line studies. The data used to derive the CIMFs are shown with open 
symbols. The predictions of our model are plotted with dots, retaining the corresponding colors from Figs. lH4l 



allowed for study of clumps with sizes > 0.1 pc in complexes 
at distances over 1 kpc. Therefore it is instructive to com- 
pare results from such studies with our model. In view of 
the complexes' remoteness, the considered regions were not 
included in the work of LAL10 and we don't have at our 
disposal their dust extinction mapping as a tool to trace 
and fit the general cloud structure (cf. Figs. IU3I left). To 
compensate this lack, we chose star-forming regions which 
clump population has been probed also by use of molecular- 
line maps. The observational C1MF was fitted directly and 
the clump mass-size diagram was used as an additional test. 



3.2.1 Mil 

La rge part of this massive star-forming region was mapped 
by iReid fc Wilsonl J2006ah in two SCUBA bands. We use 
here the clump data from their 850 /im map because of the 
richer statistics. The fitting of the observational C1MF and 
the corresponding clump-mass size diagram are plotted in 
Fig. [5] (top). A very good fit of the C1MF is achieved from 
the model case wkin4 that points to strongly gravitationally 
bound clumps. The predicted and the observational clump 
mass-size diagrams are also in agreement (Fig.[Ul right top). 
Note that we assumed a higher gas temperature T = 20 K 
which is about the average clump temperature derived by 

the authors. 

Although IReid fc Wilsonl (|2006ar > suggest a two power- 
law C1MF with M c h ~ 8 Mq, it seems that their result 
hints rather at an order-of-magnitude higher characteris- 
tic mass. This is supported by the slope in the mass range 
10 < m < 100 Mq which is similar to the IM slopes found 



from the referred molecular-line studies (cf. Table fTJ). How- 
ever, only two clumps with masses > 100 Mq in the ob- 
servational sample are far from sufficient for a plausible es- 
timate of the characteristic mass. To clarify the issue, we 
apply our met hod to a C 18 Q study of the south-western 
sector of M 17 (jStutzki fc Gustenlll990l ). As demonstrated in 
Fig. [S] bottom, the results for an intermediate-mass C1MF 
with M c h ~ 100 Mq are virtually the same, fitted from 
the same model case and free parameter values and slightly 
higher temperature. Moreover, the high-mass C1MF is also 
well fitted by the time-weighted mass distribution of gravi- 
tationally unstable clumps. Neglecting the mass bin > 2000 
Msol, which contains only 2 clumps, the slope Fhm is shal- 
lower than that of the stellar IMF. 



For independent confirmation of this result, we com- 
posed a twice larger sampl e, including clump d ata from a 
similar study of NGC 7538 JReid fc Wilsonll2005r i. As shown 
in Fig. [7] the C1MF is fitted again from the model case 
wkin4 and a similar parameter set (/3, b) for a large range of 
clump masses 7 < m < 1000 Mq, applying time- weighting 
for the more massive clumps. In comparison to the mass 
distribution from the M 17 sample, the IM slope shallows 
slightly and approaches the values from the molecular- line 
studies of nearby MC complexes (cf. Table [TJ while the HM 
slope is much steeper than that from Stuzki & Glisten's 
study and is identical with the one in Orion A (Fig. [2]) al- 
though with a smaller M c h. The agreement of the model 
with the clump data on the m — r diagram is again excellent 
(Fig. m right). 
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Figure 6. C1MF and clump mass-size diagrams of the star-forming region M 17 from dust-continuum (top) and molecular-line (bottom) 
observations. The shaded areas denote the observational mass range restrained by the lower-mass confidence limit from the model - the 
corresponding data used to derive the CIMFs are shown with open symbols in the M-R diagram. The predictions of the model are plotted 
with the same colors in the left (lines) and right (dots) panels. 
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Figure 7. Combined C1MF and clump mass-size diagram of the star-forming regions M 17 and NGC7538 from dust-continuum obser- 
vations. The designations are the same like in Fig. [6] The values of k = 0.065 and T = 30 K are retained without change from the fitting 
of the M 17 sample only. 



3.2.2 Rosette molecular cloud 

The abundant observational studies of this complex in the 
outer Galaxy indicated active star formation in the past and 
nowadays. Its clump population has been studied both on 
m olecular-line and dust cont i nuum maps; we chose the works 
of IWilliams. Blitz fc Starkl (|l995l ) and |Pi Francesco et all 
(|2010l ). respectively. The high-quality Herschel data used 
by the latter authors allow for derivation of the C1MF in 
a mass range spanning three orders of magnitude (Fig. 



[8] top left). Some incompleteness is sensible in the mass 
bins > 10 Mq and has - in our view, - two possible ex- 
planations. First, for the sake of our study we selected 
only starless clumps fro m the original sample. Second, 
|Pi Francesco et all (|2010|) applie d the Getsources algo- 
rithm ( Men'shchikov et al.l 120121 ) for clump decomposition 
which favors identification of more compact (and less mas- 
sive) objects. 

Because of the lower mass-limit of confidence, our sta- 
tistical method is not able to fit the C1MF below ~ 2 Mq. 
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The observed mass distribution for larger masses could be 
interpreted as single power-law C1MF with a slope about — 1 
(not shown). In view of the possible incompleteness of the 
high-mass clump data, we suggest rather a two power-law 
mass function, fitted from the model case wkin2 (Fig. [8] top 
left). The characteristic mass for wkin2 and for the obtained 
best-fit parameters is ~ 9 Mq (see Fig. 3 in Paper I). Adopt- 
ing this value for Mch, one gets a bit steeper intermediate- 
mass C1MF than found for other MC complexes in this Pa- 
per an d a high-mass C1M F, with a slope typical for fractal 
clouds (|Elmegreenll 19971 ). On the clump mass-size diagram, 
a shift of the model from the observed clumps is evident like 
in some other studied regions: by a factor of 2 to 4 (Fig. 
[8j top right; cf. Fig. [5]). Nicola Schneider (private commu- 
nication) provided for us estimates of the cloud mass M as 
a function of spatial scale L, calculated from the original 
Herschel column density map of Rosette. That enabled a 
check of the derived best-fit model with the observed cloud 
structure, in terms of the LAL10 work (Figs. [Tl3l left). The 
shape of the modeled curve x(L) is consistent with the dust- 
continuum data but with a shift toward lower x. Decreas- 
ing the mapping resolution parameter down to re — 0.02, 
we found excellent agreement with the data (Fig. [HI right, 
embedded diagram). That is illustrative how the observed 
general cloud structure affects the re-parameter space of the 

C1 MF fitting. 

I Williams. Blitz fc Starkl (119951 ) derived the observa- 
tional C1MF from CO mapping of Rosette MC. As seen in 
Fig. [8j (bottom), their c lump size range practicall y does not 
overlap with the one of lDi Francesco et al.l (|2010l ) while the 
mass range spans from 10 to > 1000 Mq. The C1MF is 
well fitted from the model case wkin2mag and a very 'soft' 
velocity scaling, j3 ~ 0.20. Such values of the scaling in- 
dex were found from magneto-hydrodynamic simulations of 
ICollins et all J2012t ). 



4 DISCUSSION 

The comparison with the considered molecular-line and 
dust-continuum studies of Galactic MC complexes demon- 
strates the applicability of a statistical approach to derive 
the C1MF presented in Paper I. The results on fitting the 
observational C1MF through our model are summarized in 
Table [5] They are obtained by use of two different primary 
fitting criteria: fitting of the general cloud structure x(L) — L 
or direct fitting of the C1MF. Regardless of the type of obser- 
vational data, the applied fitting criterion and of the variety 
of model cases, the best-fit parameters span relatively nar- 
row ranges: 

• 'Soft' velocity scaling: 0.20 < /3 < 0.46 

• Mainly solenoidal forcing: 0.30 ^ b ^ 0.46 

• Mapping resolution from few percent to one tenth of 
the spatial scale: 0.02 ^ re ^ 0.10 

• Typical temperatures for molecular gas phase: 8 ^ T ^ 
35 K. 

A small velocity scaling index /3, close or equal to the Kol- 
mogorov value for incompressible turbulence (0.33), may 
seem unrealistic in view of the high compressibility of in- 
terstellar turbulence which implies /3 ~ 0.50. In fact, nu- 
merical simulations of magnetized clouds show that the 



velocity power spectru m can be ev e n sha llower than in 
the Kolmogorov theory. I Collins et al.l (|2012f ) measured /3 = 
0.23 — 0.29 for thermal-to-magnetic pressure ratio in range 
0.2 — 2.0. (In our modeling, the latter value could be as low 
as ~ 0.05 in case wkin2mag; see Eq. [3]) Apparently, the ve- 
locity scaling index for incompressible turbulence /? = 0.33 
should not be necessarily treated as the lowest possible value 
from theoretical considerations. 

The best-fit range of the turbulent forcing parameter 
b is narrower, with typical values 0.38 — 0.40, in most con- 
sidered complexes. Since b can vary within a star-forming 
region, this result might be explained as a statistical effect. 
Indeed, 6 ~ 0.40 corresponds to a natural mixture between 
solenoidal and compressive modes as the latter represent 
longitudinal waves, occupying one of the three spatial di- 
mensions (see iFederrath et alj|20ld . and Fig. 8 there). The 
typical best-fit values of the mapping resolution parameter re 
are about several percent. These are the expected values, ap- 
propriate to distinguish substructures which are significantly 
smaller than the spatial scale L and significantly larger than 
the scale of dissipation. 

Generally, our results lend support to a two power-law 
shape of the C1MF: intermediate-mass and high-mass part, 
with two distinct values of the characteristic mass M c ^: ~ 10 
and ~ 200 M Q . (The molecular-line study of Rosette MC 
with its high M c h is the only exception.) The first value 
is consistent with model cases of a "virial-like" equipar- 
tition between gravitational and turbulent energy, possi- 
bly with contribution of the thermal energy at small scales 
and negligible contribution of the magnetic energy (wkin2, 
wkin2th2). These equipartition relations are not sufficient 
to argue that the modeled clumps are in virial equilibrium 
(|Ballesteros-Paredesl 120061 ) but rather indicate their gravi- 
tational boundedness or contraction. On the other hand, if 
the considered model cases hold fo r a whole cloud, they are 
indic ative for its global collapse (|Vazquez-Semadeni et al.l 
120071 ). Large characteristic masses of hundreds Mq are ob- 
tained for strongly gravitating clumps, with possible mag- 
netic support (wkin4, wkin2mag). Such are evidently the 
cases in Taurus, Orion A, M 17 and NGC 7538. 

The slope of the modeled intermediate-mass C1MF is 
shallow and does not vary significantly from complex to com- 
plex (-0.26 > r;^ d ^ -0.56) while the variety is larger for 
the high-mass one: from typical slopes for fractal clouds to 
slopes, a bit steeper than the one of the stellar IMF (—0.9 > 
Fhm > —1.5). It seems that the single power-law observa- 
tional CIMFs of slope -0.7 > F o bs > -0.8 derived from 
molecular- li ne studies in 1990s jStutzki fc Glisten! 1 19901 ; 
iBlita 1 19931 : iHeithausen et al.l 1 19981 ) have been products of 
a combination between intermediate- and a few bins from 
high-mass C1MF. On the other hand, the variety of high- 
mass C1MF slopes probably reflects a real variety of physi- 
cal conditions in individual complexes and the gravitational 
balance and evolution of the dense fragments in them. Its 
model slope V§m similar to that of the stellar IMF when the 
C1MF is time- weighte d. When derived from molecular-line 
mapping of Orion A (jTatematsu et al.lll993|) or the dust- 
continuum st udies of M 17 and NGC 7538 bv lReid fc Wilsonl 
(|2005l . l2006al ). F^m exceeds noticeably the Salpet er value 
but l ies still within the observed variability range (jKroupal 
l200ll ). 
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Figure 8. C1MF and clump mass-size diagrams of the Rosette MC from dust-continuum (top) and molecular-line (bottom) observations. 
The designations are the same like in Fig. [5] A comparison between the predictions of the illustrated models and the observed structure 
from a column density map using Hcrschel FIR-continuum data (N. Schneider, priv. comm.) is embedded in the top-right panel. 



The dust extinction mapping of Taurus, Orion A and 
Orion B by LAL10 enabled us to make a link between the 
predicted general cloud structure and the C1MF. We point 
out also the excellent agreement applying this criterion to 
the dust-continuum study of Rosette, making use of a col- 
umn density map of Rosette (Fig. [8] top right, embedded) 
obtained from Herschel data (see ISchneider et al.l |2012| . 
for details). (Note that the exti nction map derived from 
near-IR extinction using 2MASS (|Schneider et al.ll201l1 ) de- 
livers similar values for Ay < 10 m .) Generally, good agree- 
ment is found for the chosen scales of clump generation 
L between ~ 0.5 and ~ 15 pc while in Taurus the up- 
per limit is about half of that. Those are the conservative 
limits of the inertial range of turbulencqfl that is a corner 
stone in our framework of clump description. The discrep- 
ancies for larger L are to be expected and can be inter- 
preted as reflecting changes in the physical conditions - 
such scales may lie outside the inertial range and/or the 
assumption of isotherm ality does not hold for them (e.g. 
iHennebelle et al.l 120081 ). On the other hand, at L < 1 pc 
the density PDF deviates from the lognormal sh ape and 
develops a power law tail in the high-d ensity part (jKIessenl 
l200d ; lKritsuk. Norman fc WagneifeOllI ). The predicted typ- 
ical clump size in such density regimes approaches ~ 0.1 pc 
which is about the sonic scale for temperature range 10 — 
20 K, i.e. the assumption for supersonic turbulence as a gen- 
erator of clumps breaks down. 



We recall here the comment in Scct, l2~Tl i) 



Inclusion of data from dust-continuum studies in the 
analysis of the clump mass-size (771 — r) diagrams confirms 
the results from molecular-line mappings (Sect. [3~1.5|) : con- 
sistency between the slopes within 2<r and, occasionally, a 
shift by a factor of up to 3 on the size axis. Generally, the 
agreement between the locations of the 'average clump en- 
semble' and the observed clumps is better; especially, for 
the combined clump sample from M 17 and NGC 7538. The 
discrepancies on the m — r diagrams could be attributed 
to the various clump-finding techniques used in the referred 
observational studies (cf. Table [T] column 10) and hence to 
different clump size definitions. The latter would not lead 
to substantially different mass estimates if a typical clump 
structure is a dense central region and diffuse outer shell. 
A careful comparative analysis of the basic clump-finding 
algorithms, applied to identical observational or numerical 
datasets, is necessary for a more comparison between our 
statistical clump description and the observed clumps in MC 
complexes. 

The role of the magnetic fields in shaping the physical 
characteristics of MCs and of their substructures can be ac- 
count ed for more t h oroug hly in a further extension of this 
work. iMolina et al.l (|2012i ) proposed an analytical model of 
the relation between the width of the PDF a in magnetized 
turbulent medium. R includes a modification of equation [SI 
consistent with a scaling law of magnetic field B oc (p) 1 ' 2 
(equation[3]) and can be easily incorporated in our modeling. 

Finally we caution the reader that the predictions of our 
model hold for less evolved clumps of sizes > 0.1 — 0.2 pc 
and must not be compared with observational mass distri- 
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Table 2. Summary of the results on fitting the observational CIMF in the considered molecular cloud complexes. The probed parameter 
ranges (below the table) and the ranges of the model parameters that give good fits (Columns 3-7) are specified. Notation: FC = Fitting 
Criterion (x — L - general cloud structure on x(L) — L diagram; CIMF - direct fitting of the CIMF), MC = model case. 



SF region 


FC 


MC 


Model parameter ranges 
13 b k 


T[K] 


pobs 
1 IM 


CIMF 

pmod pobs 
1 IM 1 HM 


pmod 
1 HM 


m — r slope 

-,obs ^mod 


Molecular line studies 






















Taurus 


x — L 


wkin4 


0.30-0.36 


0.38-0.42 


0.07-0.09 


10-20 


-0.27±o.i7 


-0.46±0.09 


- 


- 


1.95±0.23 


1.75 


Orion A 


x — L 


wkin2mag 


0.38-0.42 


0.38-0.42 


0.07-0.09 


18-22 


-0.14±0.09 


-0.32±o.i6 


-1.82±0.13 


-1.53±o.03 


1.14±0.40 


1.42 


Orion B 


x — L 


wkin2 


0.30-0.36 


0.30-0.36 


0.055-0.065 


10-15 


-0.17±o.io 


- 


-1.28±o.l8 


-1.26±0.08 


1.95±o.i5 


1.62 




x — L 


wkin2th2 


0.30-0.36 


0.30-0.36 


0.050-0.065 


10-15 


-0.17±o.io 


- 


-1.28±0.18 


-1.16±0.05 


1.95±o.i5 


1.60 


MC sample 


CIMF 


wkin2 


0.30-0.36 


0.42-0.46 


0.02-0.10 


10-12 


— 0.44±o.io 


-0.43±o.o7 


-1.77±0.15 


-1.20±o.io 


2.09±o.i6 


1.60 




CIMF 


wkin2th2 


0.42-0.46 


0.36-0.40 


0.02-0.10 


8-10 


-0.25±o.03 


-0.32±o.n 


-1.45±0.03 


-1.24±0.14 


2.09±o.i6 


1.65 


M17 


CIMF 


wkin4 


0.43-0.45 


0.37-0.40 


0.02-0.10 


28-35 


-0.34±o.05 


-0.31±o.n 


-1.03±o.2o 


-1.27±0.02 


1.77±o.i9 


1.82 


Rosette 


CIMF 


wkin2mag 


0.18-0.23 


0.30-0.36 


0.02-0.10 


8-15 


-0.28±o.03 


-0.26±o.o9 


-1.13±0.24 


-1.06±o.io 


2.40±o.i7 


1.30 


Dust continuum studies 






















M17 


CIMF 


wkin4 


0.40-0.44 


0.37-0.40 


0.02-0.10 


28-35 


-0.35±o.05 


-0.37±o.(>6 


- 


- 


2.13±o.3() 


1.78 


M17+ 


CIMF 


wkin4 


0.43-0.46 


0.38-0.42 


0.02-0.10 


28-35 


-0.25±o.o4 


-0.33±o.io 


-1.82±o.n 


-1.25±0.09 


2.05±0.08 


1.89 


Rosette 


CIMF 

(x-L) 


wkin2 


0.44-0.46 


0.41-0.45 


0.02-0.03 


18-25 


-0.57±o.i3 


-0.56±o.o7 


-1.04±0.33 


-0.93±o.02 


2.22±o.25 


1.86 



Probed parameter spaces: 0.18 ^ /3 ^ 0.65; 0.30 < b ^ 0.55; 0.02 ^ k ^ 0.10; 8 ^ T s£ 40 K 



butions of dense (prestellar) cores. The latter ones will be 
subject of another study which consider s the high-density 
power-law tail of the density PDF (cf . iKainulainen et al.l 
12009 : iKritsuk. Norman fc WagneifeOllf ) in active star form- 
ing regions as a basis for cor es' description. The essence of 
this a pproach is presented in lDonkov. Stanchev fc Veltchevl 
(|2012j ) while the model is in process of development. Such 
study may answer the question whether the difference be- 
tween the clump mass function and the core mass function 
is physical or statistical. 



5 SUMMARY 

Our statistical approach for physical description of conden- 
sations (clumps) formed through a turbulent cascade dur- 
ing the early MC evolution predicts: i) the cloud structure 
in terms of effective size-mass scaling relations, through the 
mass-density exponent i as a function of the spatial scale 
L (DVK11); and ii) the composite clump mass function 
(CIMF; Paper I). Different models within this framework are 
generated by choosing an appropriate energy equipartition 
relation and a set of 4 free parameters: velocity scaling index 
/3, turbulent forcing parameter b, mapping resolution param- 
eter k and temperature T. In this Paper we compared CIMFs 
from molecular-line and dust-continuum studies of Galactic 
cloud complexes with ones derived from our model applying 
alternative fitting criteria: fitting the structure x(L) of indi- 
vidual complexes (when additional dust-extinction data are 
available) or direct fitting of the observational CIMF. 

Both fitting criteria lead to modeled clump mass distri- 
butions, in good or excellent agreement with the considered 
observational data. The equipartition relations which yield 
these fits indicate gravitational boundedness of the domi- 
nant clump population in the considered clouds, possibly 
including the contribution of magnetic or thermal energy: 



model cases wkin2, wkin2th2 or wkin2mag. On the other 
hand, the results for Taurus, M 17 and NGC 7538 rather 
hint at strongly gravitating or contracting inner parts of 
those complexes (model case wkin4). The derived best-fit 
values of the parameters ft, b and T for all studied individual 
clouds span relatively narrow ranges. In most MC complexes 
the typical velocity scaling index is found to be similar to 
the original "Larson's first law" (Eq. [TJ (3 = 0.38) or to 
f3 ~ 0.43 t estified from recent observations (|Padoan et al.l 
l200rj , 120091 ). The best-fit values of the turbulent forcing pa- 
rameter concentrate around b ~ 0.40 which corresponds to a 
natural mixture betwe en compressive and solenoidal modes 
JFederrath et al.ll2010T l. 

The modeled clump mass distributions support a CIMF 
which might be represented as a combination of two-power 
law functions. The latter are separated by a characteristic 
mass M c h that varies typically within one order of magni- 
tude or more: from about ten to hundreds solar masses. The 
slope of the intermediate-mass CIMF is shallow and nearly 
constant: —0.25 > Tim > —0.55. The high-mass part of the 
CIMF corresponds to gravitationally unstable clumps in all 
considered model cases and hence a more appropriate de- 
scription should take into account the dynamical clump evo- 
lution. Such description is achieved through time-weighting 
of the clump mass distribution. (The effect will be similar if 
the CIMF is derived from a statistically significant sample 
containing clumps in MC complexes at different evolution- 
ary stages of clump formation.) We obtained slopes within a 
broader range —0.9 > Thm > —1.6 that inclu des the typical 



value for fractal clouds —1 (lElmegreen 19971) as well that of 
the stellar initial mass function ( Salpeterill955l ). 



Comparison between the observational and the modeled 
clump mass-size diagrams reveals agreement of the slopes 
within the 2a limit and, in most cases, a systematic shift 
toward smaller or larger clump sizes. The latter is to be 
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attributed to variety of size definitions used in the different 
clump- finding techniques. 
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